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ABSTRACT 

The  growing  application  of  computational  aerody¬ 
namics  to  nonlinear  rotorcraft  problems  is  out¬ 
lined,  with  particular  emphasis  on  the  development 
of  new  methods  based  on  the  Euler  and  thin-layer 
Navier-Stokes  equations.  Rotor  airfoil  character¬ 
istics  can  now  be  calculated  accurately  over  a 
wide  range  of  transonic  flow  conditions.  However, 
unsteady  three-dimensional  viscous  codes  remain  in 
the  research  stage,  and  a  numerical  simulation  of 
the  complete  flow  field  about  a  helicopter  in 
forward  flight  is  i.ot  currently  feasible.  Never¬ 
theless,  impressive  progress  is  being  made  in 
preparation  for  future  supercomputers  that  will 
enable  meaningful  calculations  to  be  made  for 
arbitrary  rotorcraft  configurations. 


I.  INTRODUCTION 

The  flow  fields  of  rotorcraft  provide  a  rich  vari¬ 
ety  of  challenging  problems  in  applied  aerodynam¬ 
ics.  It  is  well  known  that  much  of  the  flow  of 
the  rotating  blades  is  nonlinear,  three- 
dimensional,  and  often  unsteady,  with  periodic 
regions  of  transonic  flow  near  the  blade  tips,  and 
with  dynamic  stall  pockets  inboard.  The  blades 
also  shed  complex  vortical  wakes,  and  serious 
aerodynamic  interactions  may  arise  between  the 
major  rotating  and  nonrotating  components. 

For  many  years,  helicopter  engineers  have 
addressed  these  difficult  problems  with  a  mixture 
of  simplified  linear  aerodynamic  theories,  wind 
tunnel  data,  and  design  charts.  At  the  same  time, 
a  small  community  of  research  scientists  has  sys¬ 
tematically  explored  the  details  of  individual 
pieces  of  the  overall  problem,  indicated  by  the 
sketches  in  Figure  (1).  References  (1)  and  (2) 
provide  an  overall  picture  of  the  practical  side 
of  helicopter  aerodynamics,  and  References  (3), 
(A),  and  (5)  summarize  many  of  the  recent  studies 
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of  these  simpler  "building  blocks."  In  recent 
years,  computational  fluid  dynamics  (CFD)  has 
begun  to  offer  new  tools  to  the  rotorcraft  commun¬ 
ity;  and  Reference  (5),  in  particular,  emphasizes 
some  pioneering  applications  of  CFD  to  transonic 
tip  flow  on  the  advancing  side  of  the  rotor  disc. 

Computational  fluid  dynamics  has  had  a  major 
impact  on  fixed-wing  aircraft  design,  but  for  a 
variety  of  reasons,  comparable  applications  to  and 
payoffs  for  rotorcraft  lag  by  a  decade  or  more. 
Nevertheless,  the  combination  of  fixed-wing  suc¬ 
cesses  and  rapid  advances  in  modern  supercomputer 
software  and  hardware  have  begun  to  attract  the 
attention  of  both  government  and  industry  manage¬ 
ment.  Consequently,  fresh  challenges  and  new 
opportunities  are  arising  for  scientists  and  engi¬ 
neers  engaged  in  rotorcraft  research  and 
development . 

Significant  applications  of  CFD  are  already  under 
way  at  several  laboratories,  universities,  and 
heicopter  companies  including  ' avesti gations 
described  at  this  conference.  This  paper  is  not 
intended  to  be  a  review  of  these  activities. 
Instead,  it  selectively  highlights  a  few  of  the 
Joint  Army/NASA  efforts  by  the  author  and  his 
colleagues  to  develop,  adapt,  and  apply  the  latest 
CFD  technology  to  rotorcraft.  The  overall,  long¬ 
term  objective  of  this  work  is  to  develop  and 
validate  advanced  CFD  codes  for  three-dimensional, 
unsteady  viscous  flews  about  arbitrary,  elastic 
rotorcraft  configurations.  It  must  be  noted  that 
we  have  no  illusions  of  building  an  all-inclusive, 
comprehensive  analysis  system  based  on  CFD. 

Rather,  the  goal  is  to  provide  a  complete  numeri¬ 
cal  simulation  of  realistic  rotorcraft  flow  fields 
with  as  few  physical  approximations  as  possible. 
This  predictive  capability  will  help  to  reduce  the 
risks  for  new  rotorcraft  designs  and  will  provide 
improved  tools  to  help  engineers  increase  perform¬ 
ance  and  efficiency  and  reduce  noise  and 
vibrations. 

In  pursuing  this  goal ,  a  twofold  approach  is  cur¬ 
rently  being  taken  at  the  Ames  Research  Center: 

1)  Mature  CFD  methods  are  being  coupled  with 
integral  wake  models  and  simplified  structural 
models  for  engineering  design  and  analysis  using 
current  or  almost-avai lable  supercomputers. 

2)  New  CFD  codes  are  being  developed  and  validated 
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for  selective  calculations  of  rotorcraft  configu¬ 
rations  using  future  supercomputers. 

By  way  of  explanation,  the  relatively  mature  com¬ 
putational  methods  used  in  the  first  approach 
include  the  transonic  small -disturbance  and  full- 
potential  flow  solvers  that  have  been  used  rou¬ 
tinely  by  the  fixed-wing  industry  for  years,  but 
which  are  only  now  becoming  useful  to  helicopter 
engineers.  Their  approximations  may  limit  their 
scope  and  accuracy,  but  they  run  fairly  easily  and 
efficiently,  and  they  can  provide  extremely  valu¬ 
able  information  to  a  skilled  user. 

In  the  second  approach,  which  is  the  subject  of 
this  paper,  the  new  codes  and  their  associated 
technology  employ  Euler  and  Navier-Stokes  flow 
solvers.  Today  these  are  mostly  being  implemented 
in  a  pure  research  mode.  That  is,  they  are  cur¬ 
rently  nursed  along  by  CFD  specialists  for  highly- 
specialized  cases,  at  an  inordinate  expense  in  CPU 
time  on  the  most  powerful  supercomputers  avail¬ 
able,  and  the  results  may  or  may  not  be  realis¬ 
tic.  Whenever  possible,  useful  and  practical 
information  is  derived  from  these  exploratory 
results,  but  the  main  purpose  is  to  build  toward 
future  capabilities.  While  we  must  refrain  from 
unrealistic  projections,  such  as  replacing  wind 
tunnels  by  computations,  or  designing  new  rotor- 
craft  automatically  by  expert  systems,  the  fixed- 
wing  legacy  clearly  shows  that  today's  CFD 
research  codes  become  tomorrow's  design  tools. 


II.  INDIVIDUAL  COMPONENTS 

The  building-block  approach  depicted  in  Figure  (1) 
provides  a  framework  for  building  sophisticated 
CFD  codes  in  manageable  steps,  with  some  early 
applications  and  practical  dividends  along  the 
way.  This  section  describes  some  of  the  recent 
CFD  developments  at  the  Ames  Research  Center  for 
these  individual  components. 

A.  Static  Airfoil  Characteristics 

The  author  and  his  colleagues  began  their  long- 
range  program  of  applying  Navier-Stokes  algorithms 
to  rotorcraft  problems  with  the  two-dimensional 
airfoil  calculations  described  in  Reference  (6). 

In  this  study,  the  aerodynamic  section  character¬ 
istics  of  several  helicopter  profiles  were  com¬ 
puted  using  the  existing  NASA-Amcs  code  AhL2l>  uJ, 
which  solves  the  Reynolds-averaged,  thin-layer 
Navier-Stokes  equations  with  the  Baldwin-Lomax 
eddy-viscosity  model  to  approximate  boundary-layer 
turbulence  [8].  A  wide  range  of  Mach  number, 
angle  of  attack,  and  Reynolds  number  was  examined 
in  Reference  (6),  with  particular  emphasis  on 
transonic  separated  flow  conditions  relevant  to 


the  advancing  blade  in  high-speed  forward 
flight.  An  extensive  study  was  made  of  the  accu¬ 
racy  and  sensitivity  of  the  results  to  the  various 
numerical  parameters  and  approximations,  and 
favorable  comparisons  were  made  with  a  large  body 
of  experimental  data. 

Figures  (2)  and  (3)  give  representative  examples 
from  Reference  (6)  for  lift,  drag  and  pitching- 
moment  coefficients.  Of  particular  interest  to 
the  helicopter  community  is  the  variation  of 
with  Mach  number.  For  example,  the  large  negative 
peak  in  the  value  of  Cm  at  M  =  0.88,  called 
"Mach  cuck,"  is  well  predicted  by  ARC2D,  as  shown 
in  Figure  (3).  Stall  and  maximum  lift  at  low  Mach 
number  were  not  addressed,  however,  and  this 
important  problem  remains  unsolved  for  high 
Reynolds-number  flows.  However,  this  is  mainly 
because  of  the  turbulence  modeling,  and  not 
because  of  computational  barriers. 

More  recently,  Holst  [9]  examined  in  detail  the 
capabilities  and  limitations  of  APC2D,  along  with 
numerous  other  transonic  airfoil  codes,  and  he 
concluded  that  codes  of  this  type  are  capable  of 
producing  airfoil  data  that  are  as  accurate  as 
wind  tunnel  measurements.  This  new  capability  is 
beginning  to  be  exploited  by  specialists  in  the 
helicopter  industry,  as  a  tool  complementary  to 
airfoil  testing. 

B.  Airfoil-Vortex  Interactions 

The  lessons  learned  from  the  preceding  experiences 
with  the  quasi-steady  version  of  ARC2D  and  from 
the  adaptation  of  the  unsteady  transonic  small- 
disturbance  code  ATRAN2  to  inviscid  airfoil-vortex 
interactions  [ 10]  led  us  to  apply  time-accurate 
Euler  and  Navier-Stokes  algorithms  [11-13]  to  this 
two-dimensional  approximation  of  blade-vortex 
interaction.  Figure  (4)  shows  the  Navier-Stokes 
results  from  Reference  (11),  using  the  so-called 
prescribed  vortex  approach  and  unsteady  solution- 
adaptive  gridding.  In  this  example,  the  structure 
of  the  vortex  was  prescribed  and  maintained 
throughout  the  calculation,  but  its  convected  path 
in  space  was  computed  as  part  of  the  solution. 

This  technique,  which  is  somewhat  analogous  to 
shock  fitting  in  the  humu.’ioal  analysis  of  tran¬ 
sonic  and  supersonic  flows,  is  designed  to  over¬ 
come  the  tendency  of  most  numerical  schemes  to 
dissipate  artificially  t^e  s'-oen  g'-^diPits  with!" 
the  vortex. 

There  was  no  significant  shock-wave/boundary-layer 
interaction  in  the  case  shown  in  Figure  (4). 
Therefore  comparisons  with  other  numerical  tech¬ 
niques  [14]  revealed  that  the  Navier-Stokes  equa¬ 
tions  were  not  necessary  to  obtain  the  correct 
time-dependent  airloads  in  this  case  -  they  just 
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required  much  more  CPU  time  to  solve.  On  the 
other  hand,  the  solution-adaptive  grid  used  in 
Reference  ill)  provided  such  fine  resolution  near 
the  moving  shock  wave  that  a  lambda  structure  was 
observed  during  the  vortex  encounter,  whereas  this 
structure  was  not  revealed  in  other  calculations. 

The  prescribed-vortex  approach  of  References  (10), 
(11),  and  (  1  li )  is  appropriate  to  cases  in  which 
the  core  of  the  vortex  bypasses  the  airfoil  or 
rotor  blade,  but  it  cannot  treat  a  direct,  head-on 
encounter  that  would  significantly  alter  the  vor¬ 
tex  structure.  However,  Rai  [13]  has  developed  a 
vortex-capturing  method  with  very  low  numerical 
dissipation  for  such  cases,  using  a  higher-order 
differencing  algorithm  for  the  Euler  or  Navier- 
Stokes  equations  and  a  special  grid  structure. 
Rai's  method  is  also  suitable  for  capturing  weak 
disturbances,  such  as  acoustic  waves. 

In  fact,  a  major  motivating  factor  in  studying 
airfoil-vortex  interactions  is  the  acoustics  of 
helicopter  blade-vortex  interaction,  and  a  central 
question  here  is  the  capability  of  "conventional" 
CFD  techniques  to  capture  the  far-fLeld  acoustic 
radiation  due  to  the  nonlinear  shock-wave 
motion.  This  issue  is  being  addressed  by  James 
Baeder,  with  the  aid  of  Rai's  method  described 
above,  and  he  has  reported  his  preliminary  find¬ 
ings  in  References  (12)  and  (15).  Baeder  is  cur¬ 
rently  comparing  the  effects  of  the  numerical 
algorithms  on  the  results  in  the  near  field  of  the 
airfoil,  where  the  reactive  components  of  the 
fluctations  are  dominant,  and  in  the  intermediate 
field,  where  the  radiative  components  begin  to 
dominate.  Although  the  surface  airloads  are  not 
significantly  different,  Figure  (6)  from  Refer¬ 
ence  (12)  shows  that  the  different  levels  of 
sophistication  in  the  various  codes  produce  dif¬ 
ferences  in  the  acoustic  pressure  throughout  the 
field.  At  present,  a  fifth-order-accurate  Euler 
code  seems  to  be  the  best  overall  compromise 
between  accuracy  and  cost  of  the  computations. 
Baeder  [12]  also  studied  the  effect  of  vortex  core 
size,  miss  distance  (including  head-on  encoun¬ 
ters),  and  Mach  number,  and  he  is  currently  exam¬ 
ining  various  airfoil  shapes  to  see  whether  BVI 
noise  can  be  alleviated  by  changing  the  airfoil 
geometry.  The  results  so  far  indicate  that  the 
Mach  number  and  miss  distance  are  the  most  impor¬ 
tant  parameters  in  determining  the  radiated  noise. 

C.  Tip  Flows  anJ  Tip-Vortex  Formation 

The  formation  and  subsequent  rollup  of  tip  vor¬ 
tices  from  blade  tips  is  an  important  fundamental 
problem  that  provides  an  ideal  opportunity  to 
extend  the  capabilities  of  the  Navier-Stokes  codes 
to  a  practical  component  of  the  total  rotor  flow 
field.  It  is  well  known  that  tip  vortices  play  a 


crucial  role  in  helicopter  aerodynamics  and  that 
their  structure  is  closely  related  to  the  blade- 
tip  geometry.  However,  past  efforts  to  alter  tip 
shapes  and  vortex  structure  have  been  largely  by 
trial  and  error,  with  mixed  results,  and  r,u  con¬ 
census  has  emerged  on  the  optimum  blade  geom¬ 
etry.  This  seems  to  be  an  area  where  computa¬ 
tional  fluid  dynamics  can  improve  our  understand¬ 
ing  of  the  basic  phenomena  and  help  design  and 
analyze  better  blade  tips.  In  addition,  we  envi¬ 
sion  that  accurately-computed  vortex  structures 
will  be  used  as  inputs  to  wake  modeling  efforts 
and/or  prescribed-vortex  methods  for  calculating 
the  complete  vortical  wake  development  for  several 
rotations  of  the  blades. 

Srinivasan,  et  al.  [16,17]  have  extended  the  pio¬ 
neering  efforts  of  Mansour  [18]  and  Kaynak,  et  al. 
[19,20]  for  viscous  three-dimensional  wings  by 
concentrating  on  the  tip  region  and  the  vortex 
formation  process,  using  the  Reynolds-averaged 
thin-layer  Navier-Stokes  equations.  The  planforms 
shown  in  Figure  (7)  have  already  been  investi¬ 
gated,  and  the  British  Experimental  Rotor  Program 
(BERP)  configuration,  Figure  (8),  is  currently 
undergoing  detailed  study. 

Figure  (9)  shows  the  computed  streamline  pattern 
for  a  swept-tip  configuration  designed  by  ONERA. 
This  example  is  nonrotating  and  is  at  a  free- 
stream  Mach  number  of  0.85.  The  swept  tip  reduces 
the  shock  wave  strength,  but  it  increases  the 
concentration  of  vorticity  and  the  peak  velocities 
in  the  tip  vortex. 

Representative  results  for  a  rectangular  wing  at 
low  speeds  are  shown  in  Figure  (10)  from  Refer¬ 
ence  (16).  An  interesting  result  of  this  series 
of  calculations  was  the  sensitivity  of  the 
detailed  flow  features  to  the  geometry  of  the  tip 
cap,  and  the  importance  of  using  different  compu¬ 
tational  grid  topologies  in  modeling  the  different 
tip-cap  shapes,  e.g.  Figure  (11).  Generally,  good 
agreement  was  obtained  with  experimental  results, 
including  the  qualitative  differences  between  the 
different  tip  caps.  However,  the  peak  suction 
levels  that  were  measured  in  the  extreme  tip 
region  were  not  obtained  in  the  calculations.  The 
reason  for  this  clear  but  relatively  minor  dis¬ 
crepancy  in  the  numerical  results  has  not  been 
identified;  it  may  be  due  to  the  turbulence  model, 
to  inadequate  local  grid  resolution,  or  to  the 
thin-layer  approximation  in  the  governing 
equations. 

Another  shortcoming  of  the  tip-vortex  calculations 
done  up  to  now  is  the  distortion  and  diffusion  of 
the  vortex  structure  downstream  of  the  biade  tip, 
due  to  the  effects  of  the  sparse  grids  in  the  far 
field  and  the  numerical  dissipation  used  to 


enhance  the  stability  of  the  code.  For  this  rea¬ 
son,  the  vorticity  contours  in  the  wake  and  the 
calculated  rollup  process  can  only  be  considered 
qualitative  at  this  time.  New,  special  grid  topo¬ 
logies  and  solution-adaptive  grid  generation 
strategies  will  be  required  to  overcome  this  limi¬ 
tation,  as  discussed  in  Section  V. 

III.  ISOLATED  ROTOR  BLADES 


figure  (13)  shows  the  instantaneous  pressure  dis¬ 
tributions  computed  with  the  same  Euler  solver  for 
a  particularly  difficult  forward  flight  case.  The 
results,  particularly  the  unsteady  shock-wave 
formation  and  its  rapid  decay  in  the  second  quad¬ 
rant,  agree  better  with  the  experiment  than  in  any 
previous  numerical  study. 

B.  Results  for  a  Lifting  Rotor  in  Hover 


The  numerical  methods  of  the  preceding  section  for 
studying  blade  tip  flows  and  tip  vortices  and  the 
experience  gained  in  working  with  them  are  now 
being  applied  to  rotating  blades.  However,  the 
extension  is  not  a  small  step,  as  the  complexities 
of  the  flow  field  and  the  computer  resources 
required  are  significantly  greater  in  these 
cases.  In  this  section,  a  brief  overview  is  given 
of  some  recent  results,  current  activities,  and 
future  directions  for  several  aspects  of  flows  on 
rotor  blades. 

Two  aspects  of  our  current  research  should  be 
pointed  out.  First,  the  codes  are  being  developed 
to  treat  hover  and  forward  flight  in  a  unified 
manner,  for  arbitrary  blade  motion.  Therefore,  as 
explained  in  References  (17),  (21),  and  (22),  we 
have  chosen  to  solve  the  unsteady  Euler  and  thin- 
layer  Navier-Stokes  equations  in  an  inertial  frame 
of  reference,  rather  than  in  blade-fixed  coordi¬ 
nates.  This  results  in  some  increases  in  CPU  time 
for  hovering  cases,  where  non-time-accurate 
convergence-acceleration  techniques  might  other¬ 
wise  be  used,  but  it  allows  the  same  general  for¬ 
mulation  and  code  to  be  used  for  both  hover  and 
forward  flight. 

Second,  as  discussed  in  Section  IV,  our  long-term 
goal  is  to  include  arbitrary  rotor-body  combina¬ 
tions  and  to  calculate  the  aerodynamic  interac¬ 
tions  between  rotating  and  non-rotating  compo¬ 
nents.  This  has  influenced  our  preliminary  choice 
of  grid  topologies  to  insure  this  capability. 

Also,  the  topology  chosen  readily  lends  itself  to 
a  straightforward  treatment  of  multiple-blade 
rotors  and  to  wake-capturing  in  hover. 

A.  Results  for  Nonlifting  Hover  and  Forward  Flight 

In  the  absence  of  lift  on  the  rotor,  the  vexing 
problems  of  the  wake  are  avoided.  Figure  (12) 
from  Reference  (21)  shows  two  Euler  solutions  for 
a  six-blads  rotor  in  hover,  at  transonic  tip 
speeds  just  below  and  just  above  the  conditions 
for  the  onset  of  the  aeroacoustic  phenomenon  known 
as  delocalization  [23].  The  effects  of  blade-to- 
blade  thickness  interference  were  studied  for  two, 
four,  and  six  blades. 


As  reported  in  Reference  (22),  C.L.  Chen  recently 
computed  the  flow  for  a  lifting  two-bladed  rotor 
with  no  wake  modeling.  Instead,  the  vortical  wake 
was  captured  approximately  by  his  Euler  code 
through  the  use  of  a  periodic  boundary  condition 
between  grid  zones  containing  the  blades.  Fig¬ 
ure  (14)  summarizes  his  solution  and  shows  tne 
importance  of  this  periodic  boundary  condition. 
After  several  revolutions  of  the  blades,  solution 
converged  with  the  several  individual  tip  vortices 
more  or  less  merged  into  one  vortical  blob  that 
induced  approximately  the  correct  downwash  near 
the  blade  tips.  Consequently,  the  agreement  with 
experimental  data  in  the  tip  region  is  quite 
good.  However,  the  accuracy  of  the  solution 
deteriorated  in  the  inboard  region,  where  the 
induced  downwash  was  too  small  and  the  resulting 
blade-element  lift  was  too  large.  Much  more  work 
remains  to  be  done  on  this  problem,  but  these 
preliminary  results  with  no  wake  modeling  are 
encouraging . 

Chen  and  McCroskey  [22]  calculated  the  flow  past 
this  hovering  rotor  for  MT  =  0.44,  0.794,  and 
0.877,  with  the  best  agreement  with  the  experimen¬ 
tal  results  occurring  at  the  two  lower  Mach  num¬ 
bers.  The  results  at  Mj  =  0.877  were  in  close 
agreement  with  the  Euler  solution  of  Agarwal  and 
Deese  [24],  which  was  obtained  using  a  simple  wake 
model.  Both  Euler  solutions  produced  a  strong 
shock  wave  that  was  significantly  farther  rearward 
than  in  the  experiment.  However,  Srinivasan  and 
McCroskey  [17]  computed  this  same  case  using  a 
Navier-Stokes  code  and  Agarwal' s  wake  model,  and 
they  found  almost  perfect  agreement  with  the 
experimental  shock  wave  position  and  strength,  as 
shown  in  Figure  (15).  It  should  be  mentioned  that 
the  three  solutions  [17,21,24]  were  obtained  with 
three  different  codes  on  three  different  grids, 
and  differences  due  to  these  purely  numerical 
factors  cannot  be  discounted.  Nevertheless,  these 
are  strong-shock  conditions  that  would  be  expected 
to  produce  significant  shock-wave/boundary-layer 
interaction,  and  it  is  our  opinion  that  these 
comparisons  substantiate  this  speculation.  Of 
course,  it  remains  puzzling  that  Strawn  and 
Caradonna  [25]  achieved  almost  as  good  agreement 
with  the  experimental  results  using  a  full- 
potential  code,  neglecting  viscous  effects 
entirely . 


Srinivasan  and  McCroskey  [17]  also  examined  sev¬ 
eral  nonrotating  configurations  that  are  sometimes 
proposed  to  simulate  rotor  blade  tips  in  conven¬ 
tional  wind  tunnels.  In  this  exercise,  the  span- 
wise  distribution  of  circulation  of  the  hovering 
rotor  described  in  the  preceding  paragraph  was 
approximately  matched  in  two  nonrotating  simula¬ 
tions:  1)  by  a  linear  va.iat.ion  in  the  spanwise 

Mach  number  distribution  with  a  constant  blade 
pitch  angle  of  4.2°,  and  2)  by  a  linear  spanwise 
twist  distribution  at  a  constant  Mach  number  of 
0.877. 

The  computations  for  this  strong  transonic  case 
are  shown  in  Figure  (16).  These  results  show  that 
the  variable  Mach-number  simulation  is  a  good 
representation  of  the  rotating  blade  in  the  tip 
region.  However,  the  more  conventional  approxima¬ 
tion  of  variable  twist  at  constant  Mach  number  is 
not  even  qualitatively  correct  anywhere  on  the 
blade,  under  these  transonic  conditions.  The 
differences  inboard  are  to  be  expected,  since  the 
flow  on  the  rotating  blade  is  completely  subsonic 
there,  but  it  is  rather  surprising  that  the  nonro¬ 
tating  solution  is  so  different  in  the  tip 
region.  On  the  other  hand,  it  should  be  noted 
that  this  overall  conclusion  is  limited  to  tran¬ 
sonic  conditions;  both  types  of  nonrotating  simu¬ 
lations  agreed  well  with  the  rotating  solution  and 
with  the  rotor  test  data  at  subsonic  tip  speeds. 

In  any  case,  this  investigation  was  a  striking 
demonstration  of  the  flexibility  and  power  of  CFD 
to  gain  physical  insight,  study  novel  ideas,  and 
examine  various  possibilities  that  might  be  diffi¬ 
cult  or  even  impossible  to  set  up  in  physical 
experiments. 

IV.  ROTOR  -  BODY  COMBINATIONS 

The  aerodynamic  interaction  between  rotating  and 
nonrotating  components  is  an  important  distin¬ 
guishing  feature  of  rotorcraft  that  produces  con¬ 
siderable  additional  complications  for  CFD  analy¬ 
ses.  This  is  because  the  inner  boundaries  of  the 
computational  grids  of  most  codes  today  conform  to 
and  move  with  the  surface  of  the  respective  compo¬ 
nents,  and  the  flow  field  around  each  object  is 
contained  within  the  interior  of  one  or  more  of 
these  grid  zones.  The  art  of  transferring  numeri¬ 
cal  data  accurately  across  grid  zone  boundaries  as 
the  solution  develops  is  by  no  means  mature  for 
fixed-wing  configurations,  and  it  is  in  its 
infancy  for  problems  with  relative  motion  between 
the  various  grid  zones  around  the  individual  com¬ 
ponents.  Therefore,  much  creative  work  remains  to 
be  done  in  this  aspect  of  CFD  for  rotorcraft 
applications . 


Fortunately,  other  engineering  problems  have 
spawned  research  activities  that  are  relevant 
here;  for  example,  rotor-stator  interactions  in 
turbomachinery  and  store  separations  of  fixed-wing 
aircraft.  Representative  two-dimensional  examples 
of  grid  topologies  and  data-transfer  strategies 
for  these  two  problems  [26,27]  are  indicated 
schematically  in  Figure  (17).  The  patched-grid 
method  on  the  left,  in  which  a  single  plane  sepa¬ 
rates  the  two  zones,  seems  intuitively  attractive 
for  many  rotorcraft  configurations.  However,  the 
overset-grid  method  on  the  right,  in  which  one 
grid  system  sweeps  through  the  interior  of 
another,  may  allow  more  flexibility  for  the  grid 
topologies  for  the  individual  components  in  some 
cases.  The  idea  of  overset  grids  is  also  attrac¬ 
tive  for  capturing  tip  vortices  that  move  relative 
to  the  blade-fixed  grid  zones,  such  as  for  the 
wakes  in  forward  flight.  Both  strategies  are 
being  studied  at  the  Ames  Research  Center  and 
elsewhere. 

The  research  of  C.L.  Chen  [21,22]  described  in 
Section  III  has  been  tailored  to  the  patchcd-grid 
concept,  particularly  regarding  the  grid  topolo¬ 
gies.  This  is  illustrated  in  Figure  (18).  Here 
the  "outer"  cylindrical  grid  zones  have  identical 
circumferential  grid  lines  at  the  interface, 
thereby  minimizing  the  complexities  at  the  moving 
boundary  between  the  rotating  and  nonrotating 
outer  zones.  "Inner"  zones  that  are  more  appro¬ 
priate  for  the  local  flow  past  the  actual  blades, 
fuselage,  wing,  etc.  can  be  imbedded  within  each 
cylindrical  block.  For  example,  a  C-mesh  topology 
for  the  viscous  region  on  the  blade  and  in  the 
near  wake  is  shown.  In  such  cases,  the  interfaces 
between  the  inner  and  outer  zones  are  likely  to  be 
more  complex;  but  as  they  are  not  in  relative 
motion,  they  can  be  treated  by  conventional  zonal 
methods  [28]. 

The  implementation  of  this  overall  strategy  has 
not  yet  been  accomplished.  The  work  described  in 
References  (21)  and  (22)  is  a  limited  first  step; 
the  Euler  equations  have  been  solved  in  the  outer 
rotating  zone  only.  However,  the  inclusion  of 
inner,  viscous-layer  and  vortical-wake  zones,  as 
well  as  nonrotating  components,  has  high  priority 
in  our  research  plans. 

V.  LIMITATIONS  AND  CRITICAL  ISSUES 

Although  the  growth  trends  of  supercomputer  tech¬ 
nology  and  overall  CFD  capability  are  well  estab¬ 
lished,  it  will  be  some  time  before  the  Euler  and 
Navier-Stokes  methodologies  outlined  in  Sections 
III  and  IV  will  offer  a  significant  competitive 
edge  over  simpler  approaches.  As  noted  in  the 
Introduction,  the  two  overriding  issues  today  are, 


FLOPS  =  no.  floating-point  arithmetic 
operations  per  unit  time 


first,  that  the  CPU  time  and  memory  are  excessive 
and,  second,  that  the  accuracy  of  the  results  is 
often  questionable.  Major  improvements  in  super¬ 
computer  hardware  and  software,  numerical  grids 
and  algorithms,  turbulence  models,  and  vortical 
wake-capturing  capabilities  are  urgently  required, 
as  briefly  discussed  below.  Other  essential  cap¬ 
abilities  which  are  not  addressed  in  this  paper 
include  improved  graphical  displays,  better  pre- 
and  post-processing  of  the  enormous  data  struc¬ 
tures,  and  aerodynamic/structural  coupling. 
Finally,  there  is  an  acute  shortage  of  both 
structural-dynamics  engineers  and  technical  man¬ 
agers  who  understand  CFD. 


A.  Computer  Speed  and  Memory  Requirements 

In  addition  to  crucial  three-dimensional  effects, 
unsteadiness  is  an  important,  complicating  aspect 
of  flows  past  rotor  blades.  This  feature  is 
shared  by  the  fixed-wing  aeroelasticity  and  turbo¬ 
machinery  communities.  These  groups  have  helped 
to  extend  the  methodologies  of  quasi-steady  aero¬ 
dynamics,  generally  a  few  years  after  they  were 
first  introduced.  However,  existing  time-accurate 
codes  tend  to  have  stability  restrictions  that 
restrict  the  time  steps  to  values  which  are  much 
smaller  than  necessary  for  accurate  resolution  of 
the  relevant  unsteady  physics  of  the  flow.  As 
discussed  in  Reference  (29),  such  restrictions 
increase  the  CPU  time  by  an  order  of  magnitude  or 
more  for  Euler  and  Navier-Stokes  calculations;  and 
therefore,  they  must  be  overcome  before  complete 
rotor  flow  fields  can  be  computed  on  a  routine 
basis. 


It  is  instructive  to  examine  the  factors  that 
determine  the  CPU  time  for  various  CFD  approaches, 
which  can  be  estimated  from  the  following  formula 
from  Reference  (29): 


cpu  r  a  x  wGT  x  Ng  x  NT/FL0PS  (1) 


where:  A 


NG 


=  "numerical  inefficiency" 
factor,  >1.0 

=  number  of  floating-point 
operations  per  grid  point  per  time 
step 

=  no.  grid  points 


N-p  =  no.  time  steps 


At 


» 


=  (no.  ref.  lengths/revolution)  x 
(no.  rev.) /At 

=  nondimn.  time  step 


The  efficiency  factor,  A,  is  introduced  to  empha¬ 
size  that  the  code  may  not  take  full  advantage  of 
the  computer  being  used.  In  practice,  it  is  a 
function  of  the  programming  efficiency,  the  degree 
of  vectorization,  the  coupling  between  the  grid 
and  the  solution  algorithm,  user  experience, 
etc.  Ideally,  its  value  should  approach  unity; 
but  especially  with  the  advent  of  supercomputers 
with  novel  architecture,  it  could  be  much  larger. 

The  number  of  arithmetic  operations  per  grid  point 
per  time  step,  Wqj,  is  a  strong  function  of  the 
numerical  method;  that  is.  of  the  flow  equations, 
the  boundary  conditions,  the  solution  algorithm, 
and  the  grid.  The  quantity  NG  represents  the 
number  of  grid  points  for  a  finite-difference 
method,  the  number  of  elements  for  a  finite- 
element  method,  or  the  number  of  panels  for  a 
panel  method.  Consequently,  WGy  NG  represents 
the  number  of  arithmetic  operations  that  must  be 
performed  at  each  iteration  or  time  step,  although 
in  some  instances  with  panel  methods,  NG  log  NG 
or  NG  is  a  more  accurate  representation  than 

mg. 

Ideally  for  rotorcraft  applications,  the  total 
number  of  time  steps,  NT,  would  simply  be  the 
number  of  time  steps  per  blade  revolution  multi¬ 
plied  by  the  number  of  revolutions  needed  to 
determine  the  aerodynamic  characteristics.  How¬ 
ever,  many  nonlinear  aerodynamics  codes  have  sta¬ 
bility  or  accuracy  limits  that  are  determined  by  a 
nondimens ional  time  step,  At  =  VAt/L.  Thus  the 
maximum  permissible  value  of  At  typically 
depends  upon  the  complexity  of  the  problem,  the 
algorithm,  the  grid,  and  the  desired  accuracy. 

Finally,  the  computing  speed,  FLOPS,  is  a  function 
of  the  computer  clock  speed  and  architecture,  the 
data  management  techniques  of  the  code,  the  memory 
requirements  (in-core  or  external  memory),  and  the 
solution  algorithm.  Thus  it  is  clear  that  many 
different  factors  determine  the  CPU  time,  and 
hence  the  cost,  of  an  aerodynamic  calculation. 

In  Reference  (29),  a  detailed  analysis  was  made  of 
the  solution  times  that  would  be  required  to  per¬ 
form  fixed-wing  flutter  analyses,  using  a  wide 
range  of  contemporary  time-accurate  methods  on 
modern  supercomputers.  This  information  is  at 
least  qualitatively  relevant  to  rotorcraft  appli¬ 
cations.  Table  1,  reproduced  from  that  investiga¬ 
tion,  shows  a  breakdown  of  the  factors  in  Equa¬ 
tion  1  for  a  wing  of  moderate  complexity  under¬ 
going  3  cycles  of  oscillatory  motion  and  25  chord- 
lengths  of  travel  per  cycle,  running  on  a  computer 
with  a  nominal  sustained  rate  of  80  million 
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Table  1.  Computational  Requirements  for  Various  3-0 
Time-Accurate  Methods 


Flow  model 

WGT 

ng 

At 

CPU, 

minutes 

Memory , 
million  words 

Notes 

Nonlinear  panel 

WGT  * 

ng  3  2  : 

106 

0.05 

60 

2.0 

a,c 

Small  disturbance 

100 

105 

0.06 

8 

0.6 

b,d 

Full  potential 

600 

105 

0.04 

23 

2.0 

d,e 

Full  potential 

and  integral  B.L. 

630 

10° 

0.02 

50 

2.0 

d,e,f 

Euler 

3000 

105 

0.01 

450 

3.0 

d,e 

Euler  and  finite 

difference  B.L. 

2000 

2  X 

105 

0.01 

600 

6.0 

d.e.g 

Thin-layer 

Navier-Stokes 

3600 

10° 

0.005 

11,000 

30 

d,e 

Reynolds-averaged 

Navier-Stokes 

4500 

2  x 

10° 

0.004 

35,000 

60 

d,e 

Notes:  a. 

b. 

c. 

d. 

e. 

f. 
9- 


A  =  2 
A  =  3 

At  for  time  accuracy 

At  for  stability  limitations 

WGT  includes  100  for  grid  generation 

5%  increase  in  Wq-j.  for  boundary  layer 

WGj  =  500  in  viscous  layer,  3000  in  inviscid  region 


floating-point  operations  per  second.  Several  of 
the  additional  assumptions  are  noted  in  the  foot¬ 
note  to  Table  1.  The  solution-time  requirements 
are  compared  graphically  in  Figure  (19).  Of 
course,  all  of  these  results  are  very  approximate, 
accurate  to  one  significant  figure  at  best. 

It  should  be  mentioned  that  a  time  step  limit  of 
At  =  C  05  was  assumed  as  a  rather  subjective 
estimate  o"  what  is  required  to  resolve  accurately 
unsteady  transonic  effects  on  a  wing,  including 
significant  shock  wave  motion.  Also,  the  esti¬ 
mates  of  the  stability  limitations  on  At  are 
very  approximate,  and  the  numbers  given  are 
intended  more  to  give  a  sense  of  the  relative 
values  of  the  various  methods,  than  the  absolute 
values.  The  important  point  is  that,  the  more 
sophisticated  the  method,  the  more  severe  is  the 
stability  restriction  on  At  for  highly  nonlinear 
problems . 

Estimates  for  rotorcraft  applications  are  more 
difficult  to  determine.  For  two  revolutions  in 
forward  flight  of  a  two-blade  rotor  with  blades  of 
aspect  ratio  10  above  a  simple  fuselage,  and  for  a 
typical  implicit  thin-layer  Navier-Stokes  code 
with  algebraic  eddy-viscosity  modeling  of  turbu¬ 
lence,  the  following  values  might  be  (optimis¬ 
tically)  appropriate: 


A  =1.5 
WGT  =  4000 

ng  =  106 

At  =  0.025 

Nt  =  5000 

Then  Equation  1  yields  CPU  =  80  hours  for  a 
100-megaflop  supercomputer,  and  approximately 
30  million  words  of  memory  would  be  required  for 
this  problem.  At  the  present  time,  everyone  has 
difficulty  achieving  even  this  level  of  perform¬ 
ance  on  simpler  problems.  These  results  illus¬ 
trate  the  magnitude  of  the  challenge.  Only  the 
largest  machines,  such  as  the  Cray  2,  have  ade¬ 
quate  main  memory,  and  its  speed  is  insufficient 
for  more  than  a  few  showcase  solutions.  It  seems 
clear  that  the  megaflop  rates  of  even  the  next 
generation  of  supercomputers  will  be  a  limiting 
factor  for  practical  CFD  computations  of  helicop¬ 
ter  aerodynamics. 

B.  Improvements  in  Grid  Generation 

An  essential  step  in  solving  rotorcraft  problems 
using  computational  aerodynamics  is  the  generation 
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of  a  suitable  grid.  Although  a  significant 
advantage  of  some  panel  methods  and  of  the  tran¬ 
sonic  small-disturbance  method  is  that  the  appro¬ 
priate  boundary  conditions  can  be  applied  on  non- 
aligned  grids,  most  of  the  more  sophisticated 
methods  rely  on  body-conforming  grids.  The  zonal 
modeling,  or  multi-block,  concept  appears  espe¬ 
cially  attractive  for  complex  configurations. 

This  is  one  of  the  most  rapidly  growing  areas  of 
CFD,  and  there  have  been  several  meetings  in 
recent  years  devoted  exclusively  to  grid  genera¬ 
tion.  The  recent  surveys  by  Thompson  [28,30]  are 
particularly  noteworthy,  and  another  major  inter¬ 
national  conference  on  grid  generation  will  be 
held  in  Miami  in  December  1988. 

For  computational  efficiency,  algebraic  and  hyper¬ 
bolic  schemes  seem  better  suited  to  the  classes  of 
unsteady  problems  that  require  regenerating  the 
grid  at  each  time  step.  High  efficiency  is  neces¬ 
sary  to  exploit  the  concept  of  adapting  the  grid 
to  some  feature  of  the  solution,  such  as  cluster¬ 
ing  grid  points  in  regions  of  large  gradients,  so 
as  to  obtain  the  maximum  accuracy  with  the  least 
number  of  grid  points.  As  noted  earlier  in  con¬ 
nection  with  Figure  (4),  dynamic  solution-adaptive 
grids  have  been  successfully  used  in  two  dimen¬ 
sions,  and  comparable  capability  will  be  essential 
for  computing  three-dimensional  vortical  wakes 
accurately.  We  believe  that  the  grid  point 
requirements,  and  hence  the  CPU  times,  can  be 


reduced  by  a  factor  of  two  or  more  by  dynamic  gr-.d 
adaption.  However,  this  will  require  a  signifi¬ 
cant  and  specialized  effort. 

C.  Turbulence  Models 

The  simulation  of  the  dynamics  of  turbulence 
remains  the  foremost  challenge  in  fluid  mechanics 
today,  and  turbulence  modeling  is  probably  the 
weakest  link  in  the  chain  of  computational  aerody¬ 
namics  technology.  For  conventional  helicopters, 
retreating-blade  stall  is  probably  the  practical 
aerodynamic  problem  most  impacted  by  turbulence 
modeling,  while  for  tilt-rotor  configurations, 
accurate  predictions  of  the  separated  flow  associ¬ 
ated  with  rotor-wing  interactions  in  hover  are 
next  to  impossible  today. 

Because  of  the  limitations  of  the  computational 
power  available  in  the  foreseeable  future,  turbu¬ 
lence  modeling  throughout  the  aeronautics  commun¬ 
ity  has  taken  the  approach  of  single-point  closure 
of  the  Reynolds-averaged  t.'avier-Stokes  equations, 
and  no  single  turbulence  model  exists  that  can  be 
applied  to  a  wide  variety  of  flows.  An  overview 
of  the  range  of  possible  turbulence  models,  rang¬ 
ing  from  essentially  no  modeling  at  all,  to  the 
hypothetical  full  simulaton  of  turbulence,  is 
given  in  Table  2.  As  no  universal  turbulence 


Table  2.  Summary  of  Turbulence  Models 


Physical 

Numerical 

Model 

Generality 

Compatibility 

Remarks 

Viscous  wedge 

Very  low 

Very  high 

Shock-B.L.  interaction 

Integral  B.L. 

Eddy  viscosity 

Low 

High 

Very  good  when  highly  tuned 

Zero  equation 

Low 

High 

Needs  more  tuning 

2-equation 

Medium 

Low  to  high 

awgt  »  20J,  it  =  ? 

Reynolds  stress  equations 

High 

Low 

3-D  separation? 

Large  eddy  simulation 

Very  high 

Low  (?) 

Guidelines  for  above  models 

Complete  simulation 

Complete 

nth  generation  supercomputers 

Are  the  more  complex  models  any  better? 


simplicity  and  generality.  Some  of  the  models  in 
this  category  lead  to  stiff  equations,  and  this 
raises  again  the  problem  of  restrictive  values 
of  At  in  Equation  1 . 

From  this  brief  overview,  it  is  clear  that  turbu¬ 
lence  modeling  will  remain  a  primary  pacing  item 
in  computational  aerodynamics  over  the  next  dec¬ 
ade,  for  both  fixed-wing  and  rotorcraft 
appl ications. 


D.  Modeling  Vortical  Flows 


model  exists,  most  contemporary  researchers  are 
focusing  their  attention  on  creating  a  catalog  of 
models  based  on  fundamental  building-block  experi¬ 
ments.  Most  of  these  models  are  being  carefully 
tested  computationally  to  determine  their  capabil¬ 
ities  and  limitations. 

As  shown  in  Table  1  and  Figure  (19),  enormous 
computer  resources  are  required  to  solve  time- 
dependent  problems  with  finite-difference  simula¬ 
tions  of  the  Reynolds-averaged  Navier-Stokes  equa¬ 
tions,  even  with  simple  turbulence  models.  Fur¬ 
thermore,  even  those  solutions  that  have  been 
published  for  steaoy  flows  have  used  grids  whose 
fine  spacing  is  limited  to  the  single  direction 
nearly  normal  to  the  body,  and  hence  fall  within 
the  spirit  of  the  thin-layer  approximation.  This 
resulting  computational  process  qualitatively 
simulates  separated  flows  and  flows  with  large- 
scale  unsteady  behaviors,  but  the  accuracy  of  such 
simulations  is  still  controversial. 

At  the  Ames  Research  Center  we  are  well  aware  of 
their  limitations  and  we  are  always  looking  for 
something  better,  but  we  have  only  used  the  zero- 
equation  eddy-viscosity  model  of  Baldwin  and  Lomax 
[8]  and  the  thin-layer  approximation  for  rotor- 
craft  problems  up  to  now.  Holst  [10]  and  Kaynak 
and  Flores  [20]  showed  a  few  examples  in  which 
improved  turbulence  models  had  some  influence  on 
the  shock  wave  location  on  airfoils  and  wings,  but 
our  transonic  airfoil  studies  [6]  did  not  indicate 
a  serious  deficiency  in  the  Baldwin-Lomax  model 
for  rotorcraft  applications  below  stall.  The  main 
concern  that  we  have  uncovered  so  far  is  in  the 
details  of  the  tip  separation  and  early  stages  of 
tip-vortex  formation.  As  discussed  in  Section  II 
in  connection  with  Figure  (10),  we  not  have  been 
able  to  reproduce  quantitatively  the  measured 
local  suction  peaks  outboard  of  Y/B  =  0.96. 

In  the  larger  arena  of  applied  aerodynamic  compu¬ 
tations  for  fixed-wing  aircraft,  however,  more 
sophisticated  models  are  beginning  to  be  used  that 
are  "tuned"  in  conjunction  with  the  numerical 
procedure  for  a  specific  class  of  flow  problems. 

As  a  result,  validations  by  means  of  experimental 
comparisons  are  mandatory,  and  confidence  in  the 
absolute  values  of  the  numerical  predictions 
remains  low  for  flows  with  significant  separation, 
including  stall,  of  course. 

In  principle,  the  more  general  models  in  Table  2 
should  cover  a  wider  range  of  flows  with  less 
"tuning,"  but  helicopter  engineers  will  probably 
remain  suspicious  of  the  idea  that  "bigger  is 
better."  Nevertheless,  in  seme  rotorcraft  prob¬ 
lems  which  involve  three-dimensional  flow  separa¬ 
tion,  the  two-equation  eddy-viscosity  models  may 
turn  out  to  be  the  best  compromise  between 


Whereas  the  treatment  of  shock  waves  in  transonic 
flow  was  a  major  focal  point  for  computational 
aerodynamics  in  the  1970's,  compressible  flow 
fields  with  embedded  regions  of  concentrated  vor- 
ticity  are  gaining  prominence  today.  Again,  the 
key  role  is  being  played  by  the  fixed-wing  commun¬ 
ity.  Some  interest  comes  from  adverse  effects  of 
trailing  wing-tip  vortices,  but  the  vortical  flows 
shed  from  sharp  leading  edges  of  tactical  aircraft 
at  high  angles  of  attack  are  the  principal  driver. 

The  helicoidal  vortical  wakes  of  rotor  blades  have 
a  much  larger  influence  in  hover  and  at  low  for¬ 
ward  speeds  than  do  the  trailing  vortices  of 
fixed-wi.ig  aircraft.  These  wakes  are  complex  in 
geometry  and  structure,  and  treating  them  accu¬ 
rately  and  efficiently  is  one  of  the  greatest 
challenge-;  in  helicopter  aerodynamics  today. 
Possible  special  treatments  include  1)  coupling 
some  form  of  wake  modeling  with  the  finite- 
difference  computations,  as  in  Reference  (5), 

2)  three-dimensional  extensions  of  the  prescribed- 
vortex  or  split-potential  methods  discussed  in 
Section  II,  as  in  References  (25),  (31),  and  (32), 

3)  adapting  a  refined  computational  grid  to  the 
concentrated  vortical  regions  as  they  are  being 
computed,  4)  developing  new  vortex-preserving 
schemes  that  reduce  the  inherent  numerical  dissi¬ 
pation  in  current  codes  which  rely  on  vortex  cap¬ 
turing,  as  in  Reference  (13),  or  5)  combinations 
of  these. 


The  flow  fields  of  rotors  are  inherently  vortical 
and  three-dimensional,  and  three-dimensional  com¬ 
putations  are  inherently  expensive.  Furthermore, 
the  results  of  the  numerical  simulations  will 
generally  be  no  better  than  the  predictions  of  the 
vortical  wake  structure.  Therefore,  it  is 
extremely  important  to  establish  what  minimum 
level  of  complexity  in  the  governing  equations 
will  suffice  and  to  determine  the  most  expeditious 
way  either  to  model  or  to  capture  vortices  within 
the  computational  domains.  This  topic  is  one  of 
the  most  fruitful  and  opportune  areas  for  creative 
research,  and  it  seems  destined  to  be  an  extremely 
active  one  for  the  next  few  years. 


VI.  SUMMARY  AND  CONCLUSIONS 

The  rotorcraft  industry  appears  to  be  entering  a 
new  era  in  which  computational  fluid  dynamics  will 
play  an  increasingly  important  role  in  the  design 
and  analysis  of  advanced  aircraft.  A  major  effort 
is  unaer  way  at  the  Ames  Research  Center  to  hasten 
this  process  and  to  capitalize  on  the  arrival  of 
future  generations  of  supercomputers. 

Airfoil  codes  are  now  available  that  can  be  used 
to  complement  wind  tunnel  testing.  The  accuracy 
of  the  numerical  simulations  is  comparable  to  that 
of  airfoil  data,  but  the  entire  flow  field  struc¬ 
ture  can  now  be  readily  examined.  The  cost  for 
the  calculation  of  each  condition  is  probably 
greater  than  that  of  an  individual  test  point,  but 
the  geometry  of  the  airfoil  can  be  changed  much 
more  easily  in  the  computer  than  in  the  wind 
tunnel . 

The  major  focus  at  present  is  on  obtaining  more 
accurate  and  less  costly  three-dimensional  solu¬ 
tions.  Although  impressive  progress  has  been 
made,  a  finite-difference  simulation  of  the  com¬ 
plete  flow  field  about  a  helicopter  in  forward 
flight  is  not  currently  feasible.  The  principal 
CFD  limitations  are  the  computer  speeds  and  memory 
capacities,  algorithm  and  solution  methods,  grid 
generation,  turbulence  models,  and  vortex  model¬ 
ing.  Other  important  limitations  are  the  inade¬ 
quate  structural  and  aerodynamic  coupling,  and  a 
shortage  of  engineers  and  scientists  who  are 
skilled  in  both  CFD  and  helicopter  aerodynamics 
and  dynamics.  Nevertheless,  the  potential  bene¬ 
fits  of  CFD  to  the  rotorcraft  industry  are  even 
larger  than  tho  that  have  already  accrued  to  the 
fixed-wing  industry,  and  the  ever-increasing 
applications  of  this  new  tool  must  continue  to  be 
encouraged . 
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Figure  1.  Development  of  rotor  blade  aerodynamics  from  simpler  cases 
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Figure  2.  Calculated  NACA  0012  airfoil  characteristics,  Re  =  6  «  10^ 
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Figure  9.  Streamlines  and  tip  vortex  of  the  ONERA  swept-tip  blade  at  M  =  0.85  and  a  =  5 
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Figure  10.  Surface  pressure  distributions  and  streamline  patterns  on  an  NACA  0015  wing  at 
M  =  0.17  and  a  =  12° 
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Figure  11.  Schematic  of  tip  caps  and  grid  topologies  for  tip-vortex  studies 
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